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Abstract. Due to the lack of test-beams in ground-based 7-ray astronomy, detector 
calibration has been a major challenge in this field. However, with the use of Cherenkov 
ring-images due to cosmic-ray muons and of strong 7-ray signals, the CAT telescope 
00 \ could be rather well monitored and understood. Here we present a few oustanding 

aspects of this work. 
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00 : INTRODUCTION: THE CAT DETECTOR 

in 

O ' The CAT (Cherenkov Array at Themis) telescope records Cherenkov flashes due 

. to VHE atmospheric showers through its 17.8m^ mirror. Its camera is located 6m 

Qv^ ! from the mirror and has a 4.8° full field of view, consisting of a central region of 

546 0.12° angular diameter phototubes arranged in a hexagonal matrix and of 54 
CL|' surrounding tubes in two "guard rings" (Fig. la). Fast electronics allows a rela- 

Q . tively low 7-ray detection threshold energy of 250 GeV (at Zenith), and the fine 

grain of the camera permits an accurate image analysis. The experiment and the 
analysis method are fully described elsewhere [1,5]. Briefly, after selecting the most 
significant triggers (total charge Qtot>30 photo-electrons), good discrimination be- 
- tween 7 and hadron-induced showers is achieved by looking at the shape and the 

r> \ orientation of the images (see the events on Fig. la): since 7-ray images are rather 

c^ I thin and ellipsoidal while hadronic images are more irregular, a first cut is applied 

which selects images with a "7-like" shape; it is based on a x^ fit to a mean light 
distribution predicted from electromagnetic showers, and a probability P(x^)>0.35 
is required. Then, since 7-ray images are expected to point towards the source po- 
sition in the focal plane whereas cosmic-ray directions are isotropic, a second cut 
q;<6° is used in the case of a point-like source, where the pointing angle a is de- 
fined as the angle at the image barycentre between the actual source position in 
the focal plane and that of the image which is reconstructed by the fit ^. As a 
result, this procedure rejects 99.5% of hadronic events while selecting 40% of 7-ray 
events. Fig lb is an example of the pointing angle distribution obtained on the 
Crab nebula: the signal is clearly seen in the first bins, while a second signal can be 



The resolution per event is of the order of the pixel size, i.e. ^0.1°. 






Xotcil time OM souirce ■- 
Signol = 5814 zt 226 
Tin = 1 .72 ± 0.07' 




4 ^|Il^I|I|J|^I^I^^i Y |I^ 



ii^i^iJii^: 



Pointing A^ngle, oc ( ) 

FIGURE 1. (a) The CAT camera and three typical events: the first image is presumably due to 
a 7 shower because of its fine and cometary shape, and because it is pointing towards the source 
position at the center of the camera; the second image is certainly due to a hadron because it is 
more spread out and pointing elsewhere; finally, the third image is clearly hadronic and would be 
easily rejected, because it is signed by the presence of a muon ring, (b) Pointing angle distribution 
from a sample of data taken on the Crab nebula between 1996 and 1999, for zenith angles ranging 
from 21° to 35°. Cuts on total charge and shape have been applied (see text). The inset shows the 
ON— OFF distribution (OFF source runs are used to estimate the hadronic background): within 
6°, the total significance is 25.6(t for a ratio of durations Tqn/Toff = 2.4 (the OFF distribution 
has been renormalized). The corresponding significance for an equal amount of ON and OFF 
runs would be 32.2cr, i.e. 4.3cr in one hour. 

seen at a;~180°, due to 7-ray images whose direction has been mis-reconstructed 
by the fit 2. 

CALIBRATION OF THE DETECTOR 



Hardware monitoring 

Mrk 501 exhibited a remarkable series of fiares during the whole year 1997 [3], 
which have been very useful for detector calibration. As an example, Fig. 2a illus- 
trates the very good quality of the mechanical monitoring for data taken with a 
zenith angle between 0° and 44°: the signal appears right at the actual position of 
the source, thus validating the angular correction which is applied on each event to 
compensate for the unavoidable and zenith-dependent slight mis-alignment of the 
optical axis of the camera ^. 



^•' The global rise of the background distribution for large values of a corresponds to large 
hadronic images which were cut by the edge of the camera. This effect becomes fainter for 
larger zenith angles, since images form closer to the center of the camera; see [6] for illustration. 
■^) This is due to the weight of the camera, which bends the arms of the telescope, and to the 
telescope azimuth and altitude axis mis- alignment, especially at large zenith angles. 



'I ii ii' I " I II I I I , i' I' III 






Mm 



FIGURE 2. (a) Projected distribution (bin size=0.05°) of the reconstructed shower angular 
origins for all data of Mrk 501 in 1997: no background substraction has been performed, and 
events are selected by the shape cut only. The cross marks the actual position of Mrk 501. (b) 
Ratio of the charge integrated in each of the 546 small inner phototubes to that predicted by a 
model of muon rings. ^100 muon images have been used for these statistics. The dispersion of 
the ratio around unity is mainly due to the uncertainty on the gain values. 

Cosmic-ray muons falling onto the mirror yield ring-like images in which the light 
distribution can be easily predicted. The fine grain of the CAT camera allows a 
fine analysis of these images [4] . In this way, the overall conversion factor between 
ADC counts and incident Cherenkov photon number can be directly checked. This 
factor involves optical efficiencies, as well as phototube pedestals and gains. As 
an example. Fig 2b shows that the camera is correctly calibrated, except for a few 
channels which are not used in the analysis. In this study, particular attention 
has been paid to the wavelength-dependent aspect, by taking special runs using 
different UV-filters placed in front of the camera. 

Validation of the simulations 



On April 16*^, 1997, Mrk 501 reached ~8 times the level of the Crab nebula. 
The signal-to-noise ratio for this night is 2.7, corresponding to a 7-ray beam with 
only 30% contamination. The good atmospheric quality of this night allows the use 
of this beam for calibration through comparison with simulations. Fig. 3a shows 
the perfect agreement on the distributions of Hillas parameters [7] : it is shown for 
events selected both with the single orientation cut and with the complete selection 
including the x^ St. The distribution of the final number of pixels A^p'^ retained 
by the fit (see [5]) is also well reproduced. Furthermore, the good agreement ob- 
served on the fourth-brightest-pixel's charge Q^*^ validates the simulations at the 
trigger level, since the trigger condition requires four pixels above threshold. Fi- 
nally, Fig. 3b shows the x^/ndf and -P(x^) distributions ^, as well as the previously 
discussed parameters, expressed as functions of the total image charge Qtot: here 
again the agreement is very good. This allows the simulation to be used to calculate 



^) The P{x^) distribution is not flat, as would be expected if the x^ were performed using variables 
with Gaussian errors. The description of Cherenkov light fluctuations in showers development is 
a very difficult task (see [5] for discussion) . 
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FIGURE 3. Comparison of 7-ray shower simulations (full lines or open squares) with the data 
(ON-OFF) from the flare of Mrk 501 on April 16'^, 1997 (filled points). Statistical errors are 
negligible for simulations, (a) Distribution of several variables, within the orientation cut a<6° 
(left column) and with the cut on shape P(x^)>0.35 also (right column): Hillas parameters width 
{aw) and length (ctl), number of pixels in the image iVP™, and fourth-brightest-pixel's charge 
Q4th. ^jjj -^2^j^jjr g,j^^ P{x^) distributions, and previous parameters as functions of Logio(Qtot)- 



the 7-ray effective detection area within the selection cuts, as well as the energy- 
resolution function ^ . This is shown in Fig. 4. In particular, a clear positive bias in 
the energy reconstruction is visible for low values of the injected energy (Fig. 4b): 
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FIGURE 4. (a) 7-ray effective detection area (in lO^'m^), including the effect of event-selection 
efficiency, as a function of the energy. Each point corresponds to simulations, while full lines come 
from an analytical 2D-interpolation over energy and zenith angle; (b) Energy-resolution functions 
"^{E ^r E, cosO) vs \og{E/E), for a fixed zenith angle 9 = 30° (y-axis units are arbitrary). Each 
plot is defined for a fixed value of the real energy E, ranging from 200 GeV to 25 TeV, and runs 
on the estimated energy E. Full lines come from an analytical 3D-interpolation over the three 
variables E, E and 9. 
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FIGURE 5. 1996 to 1999 nightly integral flux of the Crab nebula above 400 GeV, from a sample 
of data with zenith angles up to 35°. The dotted lines represents the mean flux for each period, 
and the corresponding residuals are shown on the right panel. 

it is due to a trigger effect which selects those events which benefited from a pos- 
itive fluctuation of Cherenkov light during the shower development. This energy 
over-estimation disappears when going towards higher energies, where \E' becomes 
more Gaussian, with a zero mean value and a width cr~20%. 

CONCLUSION 

More details concerning calibration and spectrum measurement will be given 
in a forthcoming paper. The good stability obtained to date on the signal from 
the Crab nebula (Fig. 5), using the acceptances and energy- resolution function 
discussed above, illustrates the good quality of running conditions and the good 
understanding of the detector. Future effort will be devoted to the analysis of 
large zenith angle data (see [6] for a first study) and to the cross-cahbration with 
the Celeste experiment, operating on the same site with an energy threshold of 
~50 GeV; the recent observation of the first common events between both detec- 
tors [2] is an encouraging result in this direction. 
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